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regularity assumptions, any given ODE always admits solvable structures even though 
finding them in general could be a very difficult task. In practice a noteworthy 
simplification may come by computing solvable structures which are adapted to some 
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1. Introduction 

In recent years a renewed interest in the classical works of Lie and Cartan led to further 
developments in the geometrical methods for the study of ordinary differential equations 
(ODEs). In particular, symmetry methods for ODEs have received an increasing 
attention in the last fifteen years. It is well known that if an ODE S admits a local 
symmetry, this can be used to reduce the order of £ by one. This procedure, usually 
referred to as symmetry reduction method, is particularly useful when £^ is a fc-order 
ODE whose local symmetries form a solvable fc-dimensional Lie algebra. In this case, 
in fact, S can be completely integrated by quadratures [121 ESI ES]- However, finding 
all local symmetries of a given ODE S is not always possible. In fact local symmetries 
(classical or higher) of a fc-order ODE in the unknown u are defined by the solutions of 
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a linear PDE depending on the derivative of u up to the order k — 1. Since the general 
solution of this PDE cannot be found unless one knows the general solution of S, one 
usually only searches for particular solutions depending on derivative of u up to order 
k — 2. Therefore, in practice, £ could not be reduced by quadratures if it does not admit 
a solvable /c-dimensional algebra of such symmetries. 

Nevertheless, one may encounter equations solvable by quadratures but with a 
lack of local symmetries. Examples of this kind, well known in recent literature 
[3, El El US, [iHl El], seem to prove that local symmetries are sometimes inadequate 
and raise the question of whether a generalization of the notion of symmetry would 
lead to a more effective reduction method. For this reason, in the last decades various 
generalizations of the classical symmetry reduction method have been proposed. Among 
these, in our opinion, a special attention is deserved by the notions of X-symmetry on the 
one hand and solvable structure on the other hand. Both notions, introduced in [21] and 
^, respectively, are suitable for a wide range of apphcations [H[5l[9| [m [131 [T5|[20l [251126] 
and will be the starting point of our forthcoming discussion. In particular, we will 
provide a more effective reduction method which interconnects both notions. 

The relevance of A-symmetries is due to the fact that many equations, not possessing 
Lie point symmetries, admit A-symmetries and these can be used to reduce the order 
exactly as in the case of standard symmetries. Despite their name, however, A- 
symmetries are neither Lie point nor higher symmetries. As shown in [9], A-symmetries 
of an ODE S can be interpreted as shadows of some nonlocal symmetries. In practice 
it means that, by embedding £^ in a suitable system S' determined by the function A, 
any A-symmetry of S can be recovered as a local symmetry of S'. This interpretation 
of A-symmetries has many advantages. First of all, it makes possible to give a precise 
geometric meaning to the A-invariance property; second, it allows to reinterpret A- 
symmetry reduction method as a particular case of the standard symmetry reduction. 

For what concern solvable structures, they were introduced by Basarab-Horwath in 
[3] and further investigated in [U [5l [TTJ [26]. The original purpose was the generalization 
of the classical result relating the integrability by quadratures of involutive distributions 
to the knowledge of a sufficiently large solvable symmetry algebra. The main point of 
this extension is that the vector fields assigning a solvable structure do not need to be 
symmetries, nevertheless they allow to implement the reduction procedure. However, 
finding local solvable structures for a /c-order ODE is not much easier than finding k- 
dimensional solvable algebras of local symmetries, if any. A noteworthy simplification, in 
general, may come by computing solvable structures adapted to the admitted symmetry 
algebras, if any. It follows that one certainly takes advantage of the presence of any kind 
of symmetry. In particular, one can include the nonlocal symmetries corresponding to 
A-symmetries. 

In this paper we consider solvable structures adapted to local and nonlocal 
symmetry algebras of any order (i.e., classical and higher), and we introduce the notion 
of nonlocal solvable structure. 

The paper is organized as follows. In section 2 we collect all notations and basic facts 
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we need on the geometry of differential equations. In particular we recall the definitions 
of local and nonlocal symmetries and also the interpretation of A-symmetries as shadows 
of nonlocal symmetries. In section 3, we recall the results of [3] and [1] in a form suitable 
to our further discussion on the integration of ODEs. Then, in section 4, we discuss our 
reduction scheme of ODEs via (higher order) local and nonlocal solvable structures. In 
particular, in this section, we discuss a practical method for the computation of solvable 
structures. Finally, in section 5 we collect some examples which provide insight into the 
application of the proposed reduction scheme. 

2. Preliminaries 

Troughout the paper we assume that the reader is familiar with the geometry of 
differential equations. Nevertheless, we collect here some notations and basic facts we 
need in the paper. The reader is referred to P, [TJ [131 EHl EH EH EH] for further details. 

2.1. ODEs as submanifolds of jet spaces 

Let M and E be smooth manifolds and tt : E' — > M be a g-dimensional bundle. We 
denote by nk : ^^(vr) — * M the k-order jet bundle associated to vr and by jk{s) the /c-order 
jet prolongation of a section s of vr. Since in this paper we are only concerned with the 
case dimM = 1, we assume that M and E have local coordinates x and {x,u^, ...,u'^), 
respectively. Correspondingly, the induced natural coordinates on J''{it) will be (x, m^), 
1 < a < g, z = 0, 1, ..., k, where the ufs are defined by M°(joo(s)) = rf* (m"(s)) /dx\ for 
any section s of vr. Moreover, when no confusion arises, Einstein summation convention 
over repeated indices will be used. 

The k-order jet space J'^(vr) is a manifold equipped with the smooth distribution 
C^ of tangent planes to graphs of /c-order jet prolongations jk{s). This is the contact 
(or Carian) distribution of J^ (vr). 

In this framework a fc-th order system of differential equations can be regarded as 
a submanifold £ C J^i'n:) and any solution of the system is a section of vr whose fc-order 
prolongation is an integral manifold of the restriction C'^l of the contact distribution 
to £^. In this paper we will deal only with (systems of) ordinary differential equations £ 
which are in normal form and not underdetermined. 

The natural projections Ti^^k '■ J^i'^) ~^ J^i'^): for any h > k, allow one to define 
the bundle of infinite jets J°°(7r) — ;> M as the inverse limit of the tower of projections 
M < — E< — J\tt) < — J2(7r) < — .... 

The manifold J°°{'k) is infinite dimensional with induced coordinates (x,m^), 
1 < a < g, z = 0,1,..., and the M-algebra of smooth functions on J'^{tt) is defined 
as J-'i'H') = UiC°°(J'(7r)). The set D{7t) of vector fields on J°°{7c) has the structure of a 
Lie algebra, with respect to the usual Lie bracket [101 ESI EH]- Moreover one can also 
define a contact distribution C on J°°(7r), generated by the total derivative operator 

Dx = d^ + u^^^du-. 
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Notice that C is the annihilator space of the contact forms O^ = dul — u^j^^dx. 

Given a fc-th order differential equation £ = {F = 0}, the l-th prolongation of S is 
the submanifold S^''^ := {D^{F) = : s = 0, 1, ..., /}. Analogously, we define the infinite 
prolongation £^°°\ 

2.2. Local symmetries 

The infinitesimal symmetries of C'^ are called Lie symmetries of J^{tt)- These symmetries 
can be divided in two classes (see [23], [29]): Lie point symmetries, which are obtained 
by prolonging vector fields X on E, and Lie contact symmetries, which are obtained by 
prolonging symmetries X on J^(7r). 

By definition, symmetries of C^ which are tangent to S are called classical 
symmetries of S. 

An analogous geometric picture holds on the infinite jet spaces. However, contrary 
to the case of finite order jet spaces, symmetries of C cannot always be recovered by 
infinite prolongations of Lie symmetries. In fact it can be proved that X = ^,0^ + fJidu^^ 
is an infinitesimal symmetry of C if and only if ^, rj"- G J-'i'ir) are arbitrary functions and 

Hence, X is the infinite prolongation of a Lie point (or contact) symmetry iff ^, tJq are 
functions on E (or J^(vr), respectively). 

On J°°(7r), symmetries X of C which are tangent to £^°°^ are called higher 
symmetries of S and are determined by the condition X{F)\^(oo) = 0. 

When working on infinite jets spaces, since D^ is a trivial symmetry of C, it is 
convenient to gauge out from higher symmetries the terms proportional to D^. This 
leads to consider only symmetries in the so called evolutive form, i.e., symmetries of 
the form X = Dl.{if°')du'^, if"' := t]"- — u^^. The functions ip"" are called the generating 
functions (or characteristics) of X. 

For the applications we are interested in, we only need to consider symmetries 
of C which are tangent to 8^°^'. This choice turns out to be also convenient since it 
noteworthy simplifies computations (see [29] and [2] fore more details and other aspects 
of cxo-jets theory). 

As already remarked in the introduction, local symmetries (Lie or higher) of a 
/c-order ODE 

£ = {Uk = f{x, U, Ui, ..., Uk-l)} 

in the unknown u are defined by the solutions of a linear PDE depending on the 
derivative of u up to the order k — 1. This linear PDE is usually called the 
determining equation and is obtained by writing in coordinates the symmetry condition 
X{uk — /) 1^(00) = 0. Unfortunately, the general solution to this PDE cannot be found 
and one usually searches only for particular solutions, if any, which depends on derivative 
of u up to order k — 2. Hence, it is not unusual to run into ODEs for which we are 
unable to compute local symmetries. 
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2.3. Nonlocal symmetries and X-symmetries 

In recent literature, to overcome the difficulties due to a lack of local symmetries, various 
attempts have been made to find a more general notion of symmetry. Hence, some new 
classes of symmetries have been introduced. Among these generalizations there are 
those introduced by Muriel and Romero in [21] and known as A-symmetries. 

The notion of nonlocal symmetry we use in this paper furnishes the conceptual 
framework to deal with many of these generalizations. We follow here the approach to 
nonlocality based on the theory of coverings (see |19j and also [29j). However, since we 
only deal with ODEs, our approach will be noteworthy simplified. The interested reader 
is referred to [29] for the general theory of coverings and nonlocal symmetries. 

Definition 1 Let S be a k-order ODE on a q-dimensional bundle tt. We shall say that 
a smooth bundle k : S —^ £ is a covering for £ if the manifold £ is equipped with a 
1-dimensional distribution C = {Cp} ~ and, for any point p E £, the tangent mapping 
K^ gives an isomorphism between C and the restriction C\^^oo) of the contact distribution 
of J°^{-n) to£^°°\ 

The dimension of the bundle k is called the dimension of the covering and is denoted 
by dim(K). Below we only consider the case dim(K) = 1. 

Definition 2 Nonlocal symmetries of £ are the symmetries of the distribution C of a 

~ (oo) 

covering k : £ -^ £ 

Here follows a coordinate description of the notion of covering and nonlocal 
symmetry. Since below we will restrict our attention only to nonlocal symmetries of 
scalar ODEs, in the following formulas we will restrict only to the case when u"" = u. 

Let TT be a trivial 1-dimensional bundle over M, with standard coordinates (x,^), 
and let 

£ := {uk = f{x,u,Ui,...,Uk-i)} ■ (1) 

Below we will consider only coverings where k is a trivial bundle k : £^°°^ x M ^ £^'^\ 
hence denoting by w the standard coordinate in R, the distribution C on £^ is generated 
by the vector field 



D, ^ = D, + Hd^ (2) 

£ 

where if is a smooth function on £^°°^ x M and D^ = d^ + Uidu + ... + fdu,.^^ is the 
restriction to £^°°'^ of the total derivative operator on J°°(tc). Hence, the covering k is 
determined by the system (see [9]) 

£':={uk = f,w, = H}. (3) 

Nonlocal symmetries of £ are symmetries of the vector field ([2]) and can be 
determined through a symmetry analysis of the system (£'y°^' on J°°{jt). Therefore 
nonlocal symmetries of £ have the form 

Y = ^d, + r],du, + ij^d^^ (4) 
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with Tji = Dxirji-i) - D^{i)ui and ipi = D^{ipi_i) - WiD^{^). 

An interesting example of nonlocal symmetry occurring in literature is related to the 
notion of A-symmetry for an ODE S (see [21]). Despite their name, in fact, A-symmetries 
are neither Lie symmetries nor higher symmetries oi S. Nevertheless, as discussed in [9], 
A-symmetries can be interpreted as shadows of nonlocal symmetries. More precisely, S 
admits a A-symmetry X iS £' = {uk = f,Wi = A}, with A G C°°{S), admits a (higher) 
symmetry with generating functions of the form yj" = e^ip^^x, u, Ui, ..., Ufc-i), a = 1,2. 

Since A-symmetries are of great interest in the applications, and analogously their 
nonlocal counterparts, it is convenient to introduce the following 

Definition S If a covering system £' , defined by ^ with H = X and A G C°°{£), 
admits a nonlocal symmetry Y with generating functions 

y9" = e'"(/)o(x,M,Mi,...,Mfc_l), « = 1, 2 (5) 

then £' will be called a X-covering for S defined by (J\). 

3. Solvable structures 

In this section we will recall basic definitions and facts on (local) solvable structures in 
the form we need in our study. The reader is referred to |3j and [4] for further details. 

The original purpose of solvable structures was to generalize the classical result 
stating that the knowledge of a solvable /c-dimensional algebra of symmetries for 
an (n — A;)-dimensional involutive distribution V, on an n-dimensional manifold N, 
guarantees that V can be integrated by quadratures. Solvable structures are a 
generalization of solvable symmetry algebras which allows one to keep this classical 
result. 

Since in this paper we restrict our attention to the integration of 1-dimensional 
distributions, we will summarize the results of [3] and [Ij only in this case. 

Given a 1-dimensional distribution V =< Z >, on an n-dimensional manifold A^, 
the definition of a solvable structure for V is the following 

Definition 4 The vector fields {Yi, ..., Yn^i} on N form a solvable structure for V =< 

Z > if and only if, denoting Vq = V and V^ =< Z,Yi, ...,Yfi >, the following two 

conditions are satisfied: 

(i) P„_i = TN; 

(ii) LYj)h-i C Pft-i, V/i G {1, ...,n - 1}. 

In particular, given a solvable structure {Yi, ..., F„_i} for V, one has the flag of 
integrable distributions 

< Z >= Do C Di C ... C Vn-i = TN (6) 

with Vs-i = {[X,Y] : X,Y e PJ, and [X, Y] denoting the Lie bracket of X and Y. 

The main difference with the definition of a solvable symmetry algebra of V is that 
the fields belonging to a solvable structure do not need to be symmetries for V. This, 
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of course, gives more freedom in the choice of the fields one can use in the integration 
oiV. 

Remark 1 It is straightforward, by the definition above, that in principle a solvable 
structure for V always exists in a neighborhood of a non-singular point for V. In fact, 
if {x*} is a local chart on N such that Z = d^i, one can simply consider the solvable 
structure generated by d^j , j = 2, ...,n. This structure is in particular an Abelian algebra 
of symmetries for T>. Nevertheless, for a given distribution T>, it is difficult to find 
explicitly such a local chart. 

In the case of 1-dimensional distributions, the main result of [3] can be stated as 
follows 

Proposition! Let {Yi, ..., y„_i} be a solvable structure for a 1-dimensional 
distribution T> =< Z > on an orientable n-dimensional manifold N. Then, for any 
given volume form Q on N , V can be described as the annihilator space of the system 
of 1-forms {ui, ...,Un-i} defined as 

w. = ^(FiJ...j?.j...jF„_ijZjfi), (7) 

where the hat denotes omission of the corresponding field, j denotes the insertion 
operator (see l2^ ) and A is the function defined by 

A = Yi^Y2^...^Yn^i^ZjVt. 

Moreover, one has that 

duJn-i = 0, 

dui = mod{ui+i, ..., u;„_i} 

for any i G {1, ...,n — 2}. 

Proof. Here we only give a sketch of the proof. The reader is referred to [3] or [1] for 
more details. In view of the definition of the forms {cjj} one readily gets that 

Vs = Ann{uJs+i,---,^n-i), Vs = 0, 1, ...,n - 2. (8) 

All these distributions, for s = 0, ...,n — 1, are completely integrable for the definition 
of solvable structure. Hence, by applying Frobenius theorem to the distribution 
'Dn-2 = Ann{iJn~i) , one should get dun-i = a A Un-i for some 1-form a. On the 
other hand, by using algebraic formula for the exterior derivative d and the formula 
defining the function A, one can show that 1/A is an integrating factor for the 1- 

form Yi_i Yn^2->Z-iQ. Hence, one gets that dun-i = 0. The rest of the proof is 

based on similar computations which allow one to show that, on the integral manifolds 
of {ui+i, ...,ujn--i}, one has that 1/A is still an integrating factor for the 1-forms 
Yij...^Yij...^Yr,^ijZjn. m 

It follows that, under the assumptions of this proposition, T> can be completely 
integrated by quadratures (at least locally). In fact, in view of dun-i = 0, locally 
Un-i = dln^i for some smooth function /„_i. Now, since dujn~2 = mod{u;„_i}, uJn-2 
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is closed on the level manifolds of /„-!• Then, iterating this procedure, it is possible 
to compute all the integrals {Ji, ..., /„„i} of T) and eventually find its (local) integral 
manifolds in implicit form {/i = Ci, ...,/„_i = Cn~i}- 

Remark 2 If one knows a solvable structure, PropositionUl allows a complete reduction 
(or integration) ofD. However, ifYi is a symmetry ofD =< Z > and Z, Yi, Y2, ..., Y^-i 
is just a frame on M, then one can still define the forms {ui} andT> = Ann{uJi, ..., c<j„_i). 
In this case, one cannot completely integrate V, but just reduce its codimension by one. 
This partial reduction is a particular case of that described in the paper fTJ and in terms 
of the forms uji 's can be described as follows. Since Yi is a symmetry of V, then it is 
also a symmetry of the Pfaffian systemi generated by Ui, ...,ujn_i. It follows that locally 
the quotient of M by the action induced by Yi is a manifold M naturally equipped with 
a Pfaffian system X. Since Yijojj = 6ij, one has that X = {uj2, ...,ujn-i}- Now the above 
mentioned reduced distribution is just T> = Ann{I). 

Proposition [1] can also be applied to the integration of ODEs. In fact, in view 
of above discussion, one can think of an ODE £^ as a manifold equipped with the 1- 
dimensional distribution V =< D^ >■ In particular, equation £ defined by (^ can be 
seen as a (A; + l)-dimensional submanifold of t/^(vr) naturally equipped with the volume 
form 

Q = dx Adu A ... A duk^i, (9) 

where A denotes the usual wedge product. 

Then, by applying Proposition [1] to ODEs one readily gets the following 

Proposition 2 Let {Yi, ..., Yk} be a solvable structure for the 1-dimensional distribution 
V =< Dx > on £ defined by (QP- Then £ is integrable by quadratures and the general 
solution of £ can be obtained in implicit form by subsequently integrating the system of 
one forms ook, ■■■,^1, in the given order. 

Remark 3 The facts observed above in Remark\^ apply also to the case of ODEs. As 
an example, consider the ODE £ := {u2 = —u/x'^}. This ODE admits the classical 
symmetry xd^ ■ The reduced equation £ is equipped with the Pfaffian system I generated 
by UJ2, or equivalently by Auj2. Since Auj2 is a 1-form on the quotient £, then it can be 
written only in terms of the invariants of Yi . Now, since the invariants of Yi are just 
Ii = u,l2 = xUx, one can check that I = {dl2 — {l2 — Ii)dli/ 12} ■ Then integral manifolds 
of I are exactly the solutions of the reduced ODE dl2/dli = {I2 — Ii)/l2 obtained by 
using standard symmetry reduction method JM 



4. New applications of solvable structures to the integration of ODEs 

In this section we discuss our new reduction scheme for ODEs. 

As remarked above, solvable structures have been introduced to generalize the 
standard symmetry reduction method. Finding solvable structures, however, is not 
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always easy and in practice a noteworthy simplification may come by computing solvable 
structures adapted to an admitted symmetry algebra Q, if any. In fact, as discussed 
below, under this assumption the determining equation of a solvable structure extending 
Q become more affordable. Moreover, since the language of covering allows one to 
handle nonlocal symmetries just like local symmetries, we can take advantage of solvable 
structures adapted also to nonlocal symmetries. In this paper we consider solvable 
structures adapted to any kind (local, higher and nonlocal) of symmetry algebras. 

Now we first discuss the practical determination of solvable structures for ODEs. 
To this end, it is useful to have the following technical Lemma which expresses symmetry 
condition in the framework of multivector fields (see [ID]). 

Lemma 1 Let T> =< Xi,...,Xr > be an r -dimensional distribution on an n- 
dimensional manifold N. The vector field X is a symmetry ofD iff 

Lx{Xi)^Xl^X2^...^Xr = Q, Vi = i,2,...,r. (lo) 

In view of this Lemma, in fact, one readily gets the following 

Proposition 3 Let £ be defined by (C]), Vt be the volume form (Ej) and uJi 's defined by 
Ij^. The vector fields {Yi, ...,Yk} on S determine a solvable structure for V =< D^ > 
iff the following two conditions are satisfied: 

(a) Y,jY2j...jYkjD,M ^ 0; _ 

(b) for any s = 0, 1, ..., k — 1, defining Yq = D^, one has 

Ly,^,{Yj)AYoAYiA...AYs = 0, Vj = 0,1,...,s. (11) 

Proof, (a) By using the notations of Definition |H one has that "D^ = TS. Hence the 
vector fields Yi, ...,Yk, D^ are linearly independent and the conclusion follows, 
(b) The result follows by Proposition [1] and (fTOl) . ■ 

Remark 4 By using the volume form Vt, conditions / fli]) can also be rewritten in the 
alternative form 

Ly^^,{Y,)^ {Yo^Y,^...^Y,M) = 0, Vj = 0, 1, ..., s. (12) 

Equations ( TTTl) (or ( TT2l) ) are the determining equations for the solvable structures 
{Yi, ..., Yfc} of V =< D^ >. If one writes down all these equations in coordinate form, 
it turns out that these form a system of 

s=o V-^ + / 

nonlinear first order PDEs in the {k + l)k unknown components of the fields Yi, ..., Y^. 
It is under determined when k = 1,2 and overdetermined when A; > 3. However, this 
system has a very special form. In fact, for s = equation flTTl) gives a system 5*0 
of equations involving only the field Yi. Then, after determining a solution Yi and 
substituting in ( TTTl) . for s = 1 one gets another system Si of equations only for Y2. 
It follows that, iterating this procedure, the determination of solvable structures just 
follows by the subsequent analysis of the linear systems 5*0, ..., Sk-i- 



,fc 
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Remark 5 In view of the definition of solvable structure, any solution Yi of Sq is a 
symmetry of the vector field D^ and hence it is a symmetry of S C J^{tt). Nevertheless, 
it does not mean that it is a symmetry of the contact distribution C^ on J'^ij^)- In fact 
Yi could also be an internal symmetry of S. For the details on possible relations between 
external, internal and generalized symmetries the reader is referred to f^. 

In view of Remark [H the determining equations for solvable structures always 
admit local solutions even though they could be hardly solvable in practice. However, a 
noteworthy simplification may occur when one only knows {Yi, ..., 1^}, r < k — 1, and a 
complete system of joint invariants {71, ...,'jk-r} for the distribution Vr- In fact, using 
the same notations of Definition H] and Proposition [1] one has the following 

Proposition 4 Let V =< Z > be a 1-dimensional distribution and {Yi, ...,Yr}, 
r < k — 1, be such that LYfT>h-i C I^h-i for any h G {l,...,r}. // one knows a 
complete system of joint invariants {71, ...,7fc_r-i} for the distribution T>r, then T> is 
completely integrable by quadratures. 

Proof. Let {ji, Qj} be a local coordinate chart adapted to the leaves of the distribution 
Vr = Ann{d'ji, ...,d'jk_r-i}. It turns out that any field d^^ is a symmetry of T>r, for 
Ld^yX^lj) = 0, and hence {Fi, ..., F^, f?7i, •••, f?7j._r_i} form a solvable structure. The 
conclusion follows by a direct application of Proposition [1] ■ 

It is clear that, since the joint invariants F G C°°{S) for T>r are described by the 
overdetermined system 

D,{F)=Y,{F) = 0, Vz=l,...,r, 

the more fields Yi one knows, the more it is feasible to find a fundamental set of invariants 
{71, ..., 7fc-r}- In fact, the analysis of differential consequences of such an overdetermined 
system of PDEs often lead to the determination of the general solution. 

Above facts, in particular, apply to the case when S has an r-dimensional solvable 
algebrslil Q =< Yi, ..., Fr >m of local classical symmetries. In such cases, the problem of 
finding solvable structures which extends Q is much easier and feasible if r is sufficiently 
large. A first simple implementation of this idea, even though with a different approach, 
has been already given in the paper [17]. In fact, in that paper it is shown how such 
a kind of completion together with the results of [3] may lead to the integration of a 
second order ODE whose Lie algebra of point symmetries is 1-dimensional. 

However there is no real need to limit ourselves to classical local symmetries. In fact, 
in view of above discussion, the determination of solvable structures can be noteworthy 
simplified if one preliminarily enlarges Q by adjoining as many symmetries as possible. 
Moreover, as already observed in the previous sections, one often runs into ODEs with 
a lack of local symmetries. Hence, in such cases, to overcome the occurring difficulties 
in the reduction procedure one could also take advantage of the existence of nonlocal 
symmetries. 

I We denote by < Yi, ...,Yr >Rthc M- linear span of the fields Yi, ...,Yr- 
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However, as discussed above, nonlocal symmetries of £ can only be taken into 
account by preliminarily embedding £ in an auxiliary system £' , of the form ([3]), and 
then searching for symmetries (higher or classical) of £' . As a consequence, the nonlocal 
symmetry reduction of £ comes from the local symmetry reduction of £' . By applying 
Proposition [1] above to the distribution < Dx > on £', one readily gets the following 
analogue of Proposition [2) 

Propositions Let {Yi, ..., Yfe+i} be a solvable structure for the 1- dimensional 
distribution V =< D^. > on the covering system £' of a k-order ODE £. Then the 
distribution!? is integrable by quadratures and the general solution of £' can be obtained 
in implicit form by subsequently integrating the system of one forms ujk+i, •••,i-o'i, in the 
given order. 

Moreover, since any covering system £' of the form ([3]) is naturally equipped with 
the volume form 

Q = dx Adu A ... A duk-i A dw, 

one also has the following analogue of Proposition [3j 

Proposition 6 The vector fields {Yi, ...,Yk^i} on £' determine a solvable structure for 
V =< Dx> iff the following two conditions are satisfied: 

(a) Fijr2J...Jn+ij5,jfi ^ 0; 

(b) for any s = 0, 1, ..., k, defining Yq = D^, one has 

Ly,+AY,) AYoAYiA ... A F, = 0, Vj = 0, 1, .., s. 

Unfortunately, if one considers a general covering system £', it is not true in general 
that £' inherits the local symmetries of £. Hence, if for example £ has a solvable algebra 
Q of local symmetries, by considering a generic covering £' one not only may lose the 
local symmetries belonging to Q but also acquires one more dimension. It follows that, 
in this case, one should also determine a solvable structure with one more vector field. 

In particular, when working with solvable symmetry algebras, above facts raise 
the question of weather £' can be chosen in such a way that Q is inherited and the 
dimensional growth is compensated by the presence of some new nonlocal symmetry. 
An answer to this question is provided by the following 

Proposition 7 Let £ be a k-order ODE of the form (Cp which admits an algebra Q of 
local symmetries. The algebra Q is inherited by a covering system £' , of the form ^, 
iff H is a joint scalar invariant of Q . In particular, if d^H = 0, Q can be extended to 
the Lie algebra Q = Q ® (0^)^. If in addition Q is solvable, then also Q is solvable. 

Proof. Each X G ^ is a vector field tangent to £ of the form X = D^{il))dui- The proof 
of the first part readily follows by imposing the invariance of £^'under the vector fields 
of Q. In fact, since each X G ^ is such that X{uk — /)|^(oo) = 0, one readily gets that 
X{wi — H) = —X{H) = 0. Hence, if if is a function of the joint scalar invariants of 
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Q and dyjH = 0, one gains also the new symmetry dyj. The thesis follows by observing 
that dw commute with all the vector fields oi Q. ■ 

However, when d^jH = 0, the advantage of solvable structures on S' is that one 
does not necessarily need that the symmetry algebra Q is inherited by S'. In fact, since 
any X G ^ is such that 

for some a, then one readily gets that X is a symmetry of the distribution < D^, d^j >. 
Hence, if Q =< Fi, Y2, ...,Yr >m, one could search just for solvable structures for < D^. > 
which extend {9^, Yi, Y2, ..., Yr} and possibly include nonlocal symmetries. 

In particular all above facts apply to the special case when S' is a A-covering system 
for S. In this case, if F is a nonlocal symmetry corresponding to a A-symmetry, one 
could search for solvable structures which extend for example the 2-dimensional algebra 
spanned by dy, and Y. 

5. Examples 

The following examples provide an insight into the applications of above reduction 
scheme. All the examples make use of solvable structures adapted to symmetries of the 
given ODE. In particular, since in a solvable structure there is no need to distinguish 
between symmetries and other vector fields, all these vector fields will be treated on the 
same footing. 

Example 1. Consider the ODE 

_ uljuj - 2UU2) , . 

The algebra of classical symmetries Qo of ( TT3l) is 2-dimensional and generated by the 
functions 

(fil = XUi — u, 
Lp2 = Ml. 

Therefore classical symmetries do not suffice for the complete symmetry reduction of 



Remark 6 By reducing equation l[TS\) through the symmetries X2,Xi generated by ip2 
and v^i, in the given order, one first obtains the reduced equation U^^ = {—U^^ + 
f/"^)f/|-2^[/"^[/| and subsequently the Riccati equation Vr, = l-2r]~^V'^ -t]^^V +r]~'^V . 
Here by U = Ui,^ = u we denote the basis invariants of X2 and analogously by V = ^U^ 
and rj = U the basis invariants of Xi. Above Riccati equation can be explicitly integrated 
even though it is not possible to exhibit an explicit solution of the initial equation (E 
Below we show an alternative integration of fT^) by means of a solvable structure. 
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However, one can also search for higher symmetries. In this case, one can find that 
up to first order the algebra of symmetries is Qq © Qi where Qi is a 3-dimensional algebra 
generated by the functions 

V^s = xu{ — 2uui, 

Let us denote by X, the symmetries generated by the functions ipi above. We have that 
the only nonvanishing commutators are 

[Xi,X2] = X2, [Xi,Xi] = X4, [Xi,X5] = X5, 

[X2,X3] = — X4, [X3,X5] = X5. 

Therefore, Xi,X2 and X4 span a solvable algebra which can be used to completely 

reduce ( TT3|l . This special example shows the advantage of considering also higher order 

symmetries in the reduction of ODE. 

Here we will apply the solvable algebra {Yi = X4, Y2 = X2, Ys = Xi} in the integration 

of ( TT3|) through the subsequent integration of the forms 003, U2 and ui. 

Now, by making use of ([7]) and coordinate expressions of the fields Y^ and Z = D^, one 

gets A = {ul — 2uu2)u\ and 

2 

{u'f—2uu2)u2 ' 

, , — rl^ J- 2xU2-2ui J I {xui-2u)ui , 

{2xuui-4u'^)u'^ + {-2xuf+{2u-x)ul+2uu1)u2+ul , 
"^ {ui-2uu2)u2ul ^""1' 

1 r/ _L ''l+'"'^l"2— 2m^«2 ^ 

^ u'f—2uu2 {uf—2uu2)ufu2 

{ul-2uu2)u2 2' 

The 1-form 003 is exact and one can check that u^ = dl^ with 

Ui \2UU2 — ui 

Then, by restricting uj2 and uJi on the integral manifolds S3 := {I3 = C3}, one finds 

l + m In(u2)-e3"l 

u;2 V = dx -\ 3 -dui + ^du2., 

1 + it]^ ln(u2) — cr^ui 

WiL = —^ dui + Tr^du2. 

-^li^S 2u*U2 ^ 2^2 

In this case both u)2\-^ and tuil^ are exact 1-forms. This is due to the fact that both 
{Fi, 1^2,^3} and {^2,^1,^} are solvable structures. Then, when restricted on the level 
manifold S3, uj2 and uji are the exterior derivative of the functions 



X 



h 



m 1 (ni-l)e -1 

M2 Ul 

1- 

1 1 (l-2«i+2«f)e~ 


HI 


2m2 2 u'i 
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respectively. It follows that the general solution of (1131) can be written in the implicit 
form {Ji = ci, h = C2, I3 = C3}. 

Example 2. Consider the ODE 

As shown in [2I]|, flT^ has no point symmetries but admits a A-symmetry with A = x/u"^. 
If we consider the system 

_ x^ 1 

^1 = ^ 

a nonlocal symmetry F of f lT^ which corresponds to the A-symmetry found by Muriel 
and Romero in [21] is that generated by the functions ip^ = ue^ and ip^ = — 2e"'. We will 
search for solvable structures which extend the nonabelian algebra Q =< Yi = ¥,¥2 = 

dn, >■ 

By making use of the Maple 11 routines, it is not difficult to find the most general 
solvable structure {Fi, 1^2,^3} for ^ = (^0) (recall that, in this case, ¥q is the total 
derivative operator restricted to (TT5|) ). The determining equations (fTTj) . in fact, 
now are Ly^{¥j) A Yq A Fi A F2 = 0, j = 0,1,2, and admit the general solution 
^3 = aidx + a2<9« + a^du^ + a^dy,, with 

(4m^ + Au'^ul + Axuui + x^) F + 8u {{uui + x) 02 — ^^03) 



ai 



Sm^ + u'^iSuj + 4) + Sxuui + 2x 



2 



02,03,04 arbitrary functions of {x,u,Ui,w) and F an arbitrary function of x + 
arctan {{2uui + x)/{2u'^)). 

Nevertheless, since here we are only interested in the integration of flT4|) . we just consider 
the following particular solution (with F = 2, 02 = 04 = and 03 = — l/(2n)) 

Y3 = Ox — -—Ou^. 
2u 

In this case, since both {^1,^2, ^3} and {Y'i,y3,F2} are solvable structures, one can 

check that duj^ = duj2 = 0. Hence, by proceeding like in the previous example, one finds 

1^3 = dl^, UJ2 = dl2 with 

J2 = 2 In-u — w — ln{4u'^ul + Axuui + Au'^ + x^). 

Is = arctan f " "^^^^" ) + x. 

It follows that, when restricted on the level manifold {I2 = 02, I3 = 03}, ui is an exact 

1-form and one can check that it is the exterior derivative of the function 

2e'=2M2 

C0S(X — C3)" 

Then, one gets the general solution of flT^ in the implicit form {/i = ci, /2 = C2, /3 = C3}. 
In this case, however, one can solve with respect to u and gets the solution of (fT4l) in 
the form 



-^1 ~ 7Z \2 ~^ '^^^^ ln(cos(c3 — x)) — 2e'^^xtan(c3 — x) 



M = ± cos(x — C2) \ X tan(x — C2) — ln(cos(x — C2)) 
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This solution depends on 3 arbitrary constants, but of course one of them is inessential. 
In fact, by suitably rearranging ci one can write the general solution of flT^ in the form 

M = ± cos(a; — C2) a/Ci — xtan(x — C2) — ln(cos(x — C2)), Ci, C2 G 1 

This solution coincides with that found in [21j . 

Example 3. Consider the ODE 

_ Ml X^ + X^ X 

"^ ~ {l + x)x " A{l + x)v? " 2{l + x)u' ^ ^ 

It can be shown that flT6l) has no point symmetries but admits a A-covering with 
A = x/v?. Hence, if we consider the covering system £' defined by fITB]) together 
with Wi = A, it admits the nonlocal symmetry Y generated by ip^ = ue^ and 
99^ = — 2e"'. We will search for solvable structures which extend the nonabelian algebra 
g=^Yi = Y,Y2 = d^>. 

Contrary to the previous example, in this case it is more difficult to find the most general 
solution of the determining equations ( ITT1) . Nevertheless, one can find some particular 
solutions. In this case we just consider the following particular solution 

(1 + x)(x + 2'U'Ui)^ X + x"^ + 2uui + 2xuui 
^3 = — o — 7 — ; ^ — c'u H -—T- — ■ Ow 

6Ui{X + UUi) ^^^[X + UUi) 

As in the previous example, since both {Yi,y2?^3} and {Y"!, Y3, 1^2} are solvable 
structures, one can check that duj^ = duj2 = 0. Then, by proceeding as in the previous 
examples, one gets the general solution of ( TT6l) : 

M = ±(C-2x + 21n(l + x)) f / — ^dx] , (17) 

^ VJ (C-2x + 21n(l + x))' / 

with C G M. Notice that, as in the previous example, in the general integral of the 
covering system £' there was a third inessential constant we have gauged out in (fT7|) . 
Moreover, we have absorbed one of the arbitrary constants in the indefinite integral. 



Example 4. Consider the ODE 



2u\ 



u 

It can be shown that flTSl) has no point symmetries but admits a A-covering with 
A = a;e~^/" — 1/x. Hence, if we consider the covering system E' defined by flTSl) 
together with wi = A, it admits the nonlocal symmetry Y generated by ^p^ = xe^ 
and if'^ = xe^ . We will search for solvable structures which extend the nonabelian 
algebra Q =<Yi = Y, Y2 = d^ >. 

Also in this case, it is very complicated to find the most general solution of the 
determining equations (ITTj) . Nevertheless, one can again find some particular solutions. 
In this case we just consider the following particular solution 

Ys = u^e^du,. 



U2 

U 
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Contrary to the previous examples, {Fi, F3, Y2} is not a solvable structure, hence in this 
case only ^3 is closed. By proceeding as in the previous examples, one still can find the 
general solution of the covering system £' . Then, by solving with respect to u one gets 
the general solution of (fTSl) : 

u = ^^ ^, (19) 

with C G M. Here, as in the previous examples, we have gauged out from (fT9|) all 
inessential arbitrary constants. 

Example 5. Consider the ODE 

2"^ ' Le^l- - 1) ,, - f ^ + u\ e^l- + xu^ + ^. (20) 

\ X J \ X J X'^ 

It can be shown that (!20!) has no point symmetries but admits a A-covering with 
A = xe^/". Hence, if we consider the covering system £' defined by (1201) together 
with Wi = A, it admits the nonlocal symmetry Y generated by ip^ = e^v^ /x and 
i^P' = —e^. We will search for solvable structures which extend the nonabelian algebra 
g=<Y^ = Y,Y2 = d^>. 

Also in this case, it is very complicated to find the most general solution of the 
determining equations (ITTl) . Nevertheless, one can again find some particular solutions. 
In this case we just consider the following particular solution 

Y^ = -^dui- 
x-^ 

By proceeding as in the previous examples, one still can find the general solution of the 

covering system £' . Then, by solving with respect to u one gets the general solution of 

u = ^ ^, (21) 

C + fo+xlniy-^dx 

with C G M. Here, as in the previous examples, we have gauged out from fl2T|) all 
inessential arbitrary constants. 

6. Concluding remarks 

Solvable structures for a 1-dimensional distribution V can be determined by solving 
the system of equations fITT]) . The analysis of these determining equations is in general 
more feasible if one makes use of symbolic manipulation packages like those written for 
Maple V distributions. In particular, we suggest the packages DifferentialGeometry and 
Jets by I. Anderson (and coworkers) and M. Marvan, respectively. The first package 
is recommended for various computations involved in the determination of solvable 
structures and the further integration of the forms tUj. The second package is very 
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helpful in the computation of symmetries and the analysis of differential consequences 
of a given system of equations. 

Acknowledgement 1 The authors would like to thank G. Cicogna and G. Gaeta for 
constructive discussions and suggestions. This work was partially supported by GNFM- 
INDAM through the project "Simmetrie e riduzione per PDF, principi di sovrapposizione 
e strutture nonlocali" . 

[1] I. Anderson and M. Fels 2005 Exterior Differential Systems with Symmetry Acta Appl. Math. 87 

3-3f 
[2] I. Anderson, N. Kamran and P. Olvcr f 993 fnternal, External and Generalized Symmetries Adv. 

Math. 100 53-fOO 
[3] P. Basarab-Horwath f992 Integrability by quadratures for systems of involutive vector fields 

Ukrain. Mat. Zh. 43 1330-1337; translation in Ukrain. Math. J. 43 (1992) 1236-1242 
[4] M.A. Barco and G.E. Prince 2001 Solvable symmetry structures in differential form applications 

Acta Appl. Math. 66 89-121 
[5] M.A. Barco and G.E. Prince 2001 New symmetry solution techniques for first-order non-linear 

PDEs Appl. Math. Comput. 124 169-196 
[6] G. W. Bluman and S. C. Anco 2002 Symmetry and Integration Methods for Differential Equations 

(Berlin: Springer) 
[7] G. W. Bluman and S. Kumei 1989 Symmetries and differential equations (Berlin: Springer) 
[8] G. W. Bluman and G. J. Reid 1988 New symmetries for ordinary differential equations IMA J. 

Appl. Math. 40 87-94 
[9] D. Catalano Ferraioli 2007 Nonlocal aspects of A-symmetries and ODEs reduction J. Phys. A: 
Math Theor 40 5479-5489 
[10] S. S. Chern, W. H. Chen and K. S. Lam 2000 Lectures on Differential Geometry (Singapore: World 

Scientific) 
[11] Cicogna G, Gaeta G and Morando P 2004 On the relation between standard and /x-symmetries 

for PDEs J. Phys. A.Math. Gen. 37 9467-86 
[12] Fels M E 2007 Integrating scalar ordinary differential equations with symmetry revisited Found. 

Gomput. Math. 7 417-454 
[13] G. Gaeta 1994 Nonlinear symmetries and nonlinear equations (Dordrecht: Kluwer) 
[14] G. Gaeta and P. Morando 2004 On the geometry of lambda-symmetries and PDEs reduction J. 

Phys. A:Math. Gen. 37 6955-6975 
[15] M.L. Gandarias, E. Medina and C. Muriel 2002 New symmetry reductions for some ordinary 

differential equations J. Nonlin. Math. Phys. 9 Suppl.l 47-58 
[16] K. S. Govinder and P. G. L. Leach 1997 A group-theoretic approach to a class of second-order 

ordinary differential equations not possessing Lie point symmetries J. Phys. A 30 2055-2068 
[17] T. Hartl and C. Athorne 1994 Solvable structures and hidden symmetries J. Phys. A: Math Gen 

27 3463-3471 
[18] N. Kamran 2002 Selected topics in the geometrical study of differential equations (Providence, RL 

American Mathematical Society) 
[19] L S. Krasil'shchik and A. M. Vinogradov 1989 Nonlocal Trends in the Geometry of Differential 
Equations: Symmetries, Conservation Laws, and Backlund Transformations Acta Appl. Math. 
15 161-209 
[20] P. Morando 2007 Deformation of Lie derivative and /^-symmetries J. Phys. A: Math. Theor. 40 

11547-11559 
[21] C. Muriel and J.L. Romero 2001 New method of reduction for ordinary differential equations IMA 

Journal of Applied mathematics 66 111-125 
[22] C. Muriel and J.L. Romero and P. Olver 2006 Variational C°°-symmetries and Euler-Lagrange 
equations J. Differential equations 222 164-184 



Local and nonlocal solvable structures in ODEs reduction 18 

[23] P.J. Olver 1993 Application of Lie groups to differential equations (Berlin: Springer) 

[24] P.J. Olver 1995 Equivalence, invariants, and symmetry (Cambridge: Cambridge University Press) 

[25] E. Pucci and G. Saccomandi 2002 On the reduction methods for ordinary differential equations J. 

Phys. A 35 6145-6155 
[26] J. Sherring and G. Prince 1992 Geometric aspects of reduction of order Trans. Amer. Math. Soc. 

334 433-453 
[27] M. Spivak 1999 A comprehensive introduction to differential geometry. Vol 1 (Houston, Texas: 

Publish or Perish, Inc.) 
[28] H. Stephani 1989 Differential equations. Their solution using symmetries (Cambridge: Cambridge 

University Press) 
[29] A. M. Vinogradov et al. 1999 Symmetries and conservation laws for differential equations of 

mathematical physics (Providence, RI: American Mathematical Society) 



